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To identify the contributions of the multiband nature and the anisotropy of a microscopic electronic structure to a 
macroscopic vortex lattice morphology, we develop a method based on the Eilenberger theory near H& combined with 
the first-principles band calculation to estimate the stable vortex lattice configuration. For a typical two-band supercon- 
ductor MgB 2 , successive transitions of vortex lattice orientation that have been observed recently by small angle neutron 
scattering [Das, et ah: Phys. Rev. Lett. 108 (2012) 167001] are explained by the characteristic field-dependence of two- 
band superconductivity and the competition of sixfold anisotropy between the cr- and 7r-bands. The reentrant transition 
at low temperature reflects the Fermi velocity anisotropy of the cr-band. 

KEYWORDS: vortex lattice morphology, two-band superconductivity, MgB 2 , first-principles band calculation, 
quasi-classical Eilenberger theory 



Since vortices are self-forming objects in mixed states of 
type-II superconductors under applied magnetic fields, the 
properties of the vortex states give us valuable information 
on the characteristics of each superconductor. Among them, 
the vortex lattice morphology, depending on the temperature 
T and magnetic field H, is closely related to the anisotropy 
of the superconductivity. For example, reflecting the fourfold 
symmetric structure of a superconducting gap or the Fermi 
velocity in momentum space, the vortex lattice shows a grad- 
ual transformation from a triangular vortex lattice at low fields 
to a square lattice at high fields. 1 The anomalous vortex lat- 
tice transformation in a superconductor Nb, including a sca- 
lene triangular lattice, is considered to be a result of the 
anisotropy of the Fermi surface structure. 2_4 For a reliable the- 
oretical estimation of the vortex lattice morphology, we need 
microscopic calculations such as the use of the Eilenberger 
theory, beyond the phenomenological ones of the extended 
Ginzburg-Landau (GL) theory 45 or the nonlocal London the- 
ory. 6 Furthermore, to accurately consider the anisotropy of 
the Fermi surface structure, it is preferable that the results 
of first-principles band calculation for electronic states, such 
as density functional theory (DFT), are combined in the es- 
timation of a stable vortex lattice. This theoretical approach 
was carried out in a study of the vortex lattice morphology in 
Nb. 3 The success encouraged us to further extend the vortex 
physics along this line. 

The properties of a superconductor MgB2 as a typical 
multiband superconductor also attract much attention. The su- 
perconductivity in MgB2 consists of the cr-band with a large 
superconducting gap and the 7r-band with a small gap. 7 ' 8 In 
the vortex state, there are some interesting phenomena re- 
flecting the multiband superconductivity, such as a rapid in- 
crease in the //-dependence of low-temperature specific heat 
C(H) = y{H)T at low fields, 913 and the vortex clustering of 
the so-called type 1.5 superconductivity at low fields. 14 Also, 
an interesting vortex lattice morphology has recently been re- 
ported in MgB2 for H || c. 15 ' 16 Since the crystal lattice is 
hexagonal, a triangular vortex lattice is formed in MgB2. As 
schematically shown in Fig. 1, from small angle neutron scat- 
tering (SANS) experiments, 16 the orientation of the triangular 
vortex lattice is (p - 0° (F-phase) at low fields, and <p = 30° 



(I-phase) at high fields, where <p is the angle measured from 
the a*-axis in real space (a-axis in reciprocal space) of the un- 
derlying hexagonal crystal. In the intermediate fields, the ori- 
entation rotates gradually between 0° < <p < 30° (L-phase). 
However, at low-T, the high-field I-phase reported in a previ- 
ous work 15 has been confirmed as a metastable state in a more 
recent work. 16 The true stable state at low-T and high fields 
is found to be L-phase, as shown in Fig. 1. Therefore, a stable 
vortex lattice shows the successive transition L — > I — > L — > F 
along H C 2 upon raising T from T - to the transition temper- 
ature T c . In a previous study based on phenomenological GL 
theory, 5 the transition I — > L — > F was considered as a result 
of competition in anisotropies between the cr-band and the n- 
band in two-band superconductivity. However, the reentrant 
transition L — > I at low-T has not yet been explained theoret- 
ically. We note that the original GL theory is valid only near 
T c . To discuss the vortex morphology of MgB2, we have to 
consider the 6th- and 12th-order derivative terms for gradient 
terms of the GL equation. The nonlocal correction terms of 
the 2n-th order derivative are on the order of (1 - T/T c ) n ~ l . 
Thus, it is difficult in the framework of the GL theory to dis- 
cuss the vortex morphology in the low-temperature region, 
because we need further higher-order derivative correction 
terms. On the other hand, since the Eilenberger theory can be 
applied to all T regions, we can obtain a reliable estimation 
for a vortex structure including that in a low-r region. The 
purpose of this study is to investigate the successive transition 
ofL— > I — > L — > F along H C 2, on the basis of the Eilenberger 
theory of two-band superconductors combined with DFT cal- 
culation 17 to consider the Fermi surface of MgB2. 

First, we study the sixfold anisotropy of the Fermi surface 
structure in MgB2. The electronic state of MgB2 is evaluated 
by DFT calculation 17 using local density approximation and 
the lattice constants a - b - 3.083A, and c - 3.521A. 18 We 
chose the muffin-tin radii Ryu of 2.35 and 1.67 a.u., respec- 
tively, and use a plane-wave cutoff RMi,mmK max = 9.0. The 
first Brillouin zone is divided into 120 x 120 x 90 grids. By 
the tetrahedron method 19 for determing patches of the Fermi 
surface from the grids, we estimate the Fermi wave number 
k,, Fermi velocity v, = (v, jX , v, ?y , v !?z ), and Fermi surface DOS 
Dj(ep) on the j-th patch of the Fermi surface. 



J. Phys. Soc. Jpn. 



LETTERS 



(b) I orientation \ 




1(1 20 30 40 

Temperature (K) 



Fig. 1. (Color online) (a) Schematic phase diagram of vortex lattice mor- 
phology as a function of H and T, obtained by the SANS experiment. 16 The 
phase boundary near H c 2 is extrapolated. Triangular lattices in the right panel 
show the orientation <j> of each-phase, (b) The T dependence of <f> is schemat- 
ically presented along H^i ■ 



As shown in Fig. 2(a), our calculation reproduces the well- 
known Fermi surfaces of the cr-band and the 7r-band. 20 For 
later discussions on the reasons for the vortex lattice morphol- 
ogy in MgB2, we study the anisotropy of the Fermi veloc- 
ity projected within the xy plane perpendicular to the applied 
field along the z-direction. The color on the Fermi surfaces in 
Fig. 2(a) shows the orientation angle tp Vj of the Fermi velocity 
defined as tan</) v . = v iy /v ix . Figures 2(b) and 2(c) show the 
^,,-resolved Fermi surface DOS defined as 



D m (ip v ) = (8(<p v - tp v .))„ 



(1) 



for each band (m - cr,n), where the band-resolved Fermi 
surface average is (■■•),„ = Eiem AteX - • ■ )• As shown in 
Fig. 2(b), the Fermi surfaces of the cr-band are two hexag- 
onal cylinders rather than circular ones. Most of the Fermi 
velocity on the fiat surface is oriented to 30° or equivalent 
directions. The Fermi surface of the 7r-band is also fiat, and 
most of the Fermi velocity on the fiat surface is oriented to 
0° or equivalent directions, as shown in Fig. 2(c). Therefore, 
the D m (ip v ) values of the Fermi surface are different between 
the two bands, inducing a competition between the F and I 
phases for the vortex lattice orientation, as will be discussed 
later. The ratio of the Fermi surface DOS of the 7r-band to the 
cr-band is calculated to be N n /N cr = 1.36, indicating that the 
7r-band contribution is larger. 

Next, we explain our formulation for estimating the sta- 
ble vortex lattice configuration, which is extended to two- 
band superconductors in this work. To shorten our calcula- 
tion within a reasonable computational time even when we 
consider the detailed structure of three-dimensional Fermi 
surfaces, our study of the stable vortex lattice is restricted 
near H C 2, as a first approach to this problem. We assume 
an isotropic s-wave superconducting gap in each band, and 
a clean limit. 

Quasi-classical Green's functions g(cj n ,kj, r), /(o»„,k,-,r), 
and /'(w„,k,, r) with a Matsubara frequency a>„ are defined 
at k, on each Fermi surface. They are calculated by the Eilen- 
berger equation 13 ' 21 

{u„ + v ; • (V + iA)} f{u„, k,, r) = AJr)g(a) n , k„ r) (2) 

and the equivalent equation for /', where g = (1 — ff^) 1 ^ 2 
and A is the vector potential. Throughout this letter, we use 
Eilenberger units, 22 i.e., energy, temperature, and length are 
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Fig. 2. (Color online) (a) Fermi surfaces of MgB2 obtained in our calcu- 
lation. The two cylinders in the center are the <T-band's Fermi surface, and 
the outside rings are the ff-band's Fermi surfaces. The color on the surface 
indicates the orientation ip yj (mod 30°) of the Fermi velocity, (b) <p„-resolved 
Fermi surface DOS Do-itpy) for the cr-band. (c) D n (ip v ) for the ff-band. In (b) 
and (c), we plot lines of (D„,(p v ) cos ^..D,,,^,,) sin (/;,,) for < ip v < 360° 
(m = cr, n), so that the distances from the center to the lines indicate D m {tp v ). 



in units of nk B T c , T c , and £o = hv F /27rk B T c , respectively. The 
self-consistent gap equation for the pair potential A m (r) of m- 
bands (m - cr, n) is given by 

/ Mr) \_^ T v/ </("„, k r -,r)>, \ m 

[ A„(r) j lu\ (f(co n ,ki,r)) n J' w 



As for the interaction constants of superconductivity, 



V 



V,r 



V a 



(4) 



where Va-.a- (V„^) is the pairing interaction in the cr- (n-) band, 
and Va-,n is the Cooper pair transfer between two bands. We 
have to select these parameters under the condition to re- 
produce the gap ratio A^/A^ - \ ~ \ observed in experi- 
ments. 8 Previous first-principles band calculation 23 estimated 
(i) y<r,n = 0.15Vo. ]O - and V K ,„ - 0.32V <T - jO -. These are used in 
our calculations. We determine Va-^ from the gap equation at 
T = T c . 

We expand A m (r) and /(«„, k,, r) in terms of the Landau- 
Bloch function (/'^(r), 24,25 whose the A^-th coefficients of the 
Landau level are A m _jv and /at(w„, k,). The vortex lattice con- 
figuration is determined by the parameters of ^jv,q(r), an d tne 
relative angle of the crystal coordinate and vortex coordinate. 
Along H C 2 where g — > sgnw„, the Eilenberger equation (2) is 
reduced to 24 ' 25 



/ W (w„,k, em ) = 2_j(M l ) N ,M^rr 
M 



(5) 



with M n ,m = \u n \6 N , M + y/N + \b*5 NM -\ - y/Nbd NyM +\ and 
b - (vi X + iv;j,)sgn(to,;) V/7/(2 V2). The Fermi velocities v,-, x 
and v; v at k,, appearing in b, determine the k, dependence of 
/n. Therefore, at B — H C 2, the gap equation (3) is reduced to 



Ao-.jv 



_ p yi y / (( M 1 )n,m)v\t,m 



(6) 



In the Eilenberger theory, the free energy in a two-band 
superconductor is given by F — k 2 B 2 + F c \ + F C 2 with the GL 
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parameter k and 



F cl = J] A;(r)(y-%„A„(r), (7) 

n,m=(r,7i 

f c2 = -tY j Yj< ! ^» (8) 



where (• • • ) indicates the spatial average and I - (g - 
sgn W „)[2w„ + v,- • (Vln(///') + 2iA)] + /A* + A/'". Using 
the gap equation (3), we obtain 



Fcl = T YjYj ( A '»</ 1 ">» + A m</>'»>- 



(9) 



Since the condensation energy part F c \ + F C 2 is a simple sum- 
mation of each band contribution, the derivation of free en- 
ergy near H C 2 can be performed as in a single band case, 26-28 
as follows. We substitute g ~ 1 - ff/2 - (// t ) 2 /8 in F and 
set A = A C 2 + Ai, where V x A C 2 = H C 2, V x Aj = hi, and 
hi = H - H C 2 + h 5 for an applied field H near H Q 2- Thus, the 
Abrikosov identity 27 is derived as 

T Z Z V/ 1 ^-/ 1 + K,f)) m - 4/hiTh: = 0, (10) 

o>„>0 m 

where h, is the magnetic field induced by supercurrent, and 

w,,>0 m 



Vxh,= 



Finally, we obtain the free energy at magnetic induction B 
H + h s near H C 2 as 26 

>2 



?-*-*&?■ 



with 



r -- t-L y 2 <// t (A B /t + A;/)> m - fe. (i3) 



4k 2 M *-* n 

s a»„>U m 



m 



Since we consider the case k » 1 here, the last term with /z 2 in 
Eq. (13) is neglected. The free-energy minimum corresponds 
to the minimum of T' . 

Using Dj(ep) and v, obtained by the DFT calculation above, 
we solve the gap equation (6) of matrix form and obtain H C 2 
and A m at a given T. Substituting A m to Eq. (5) we obtain 
quasi-classical Green's function /. By these A m and /, we 
calculate T' in Eq. (13). These calculations are performed for 
possible orientations (0° < fy < 30°) of the triangular vortex 
lattice, and we determine the stable orientation of the vortex 
lattice along the H C 2 line, by the free-energy minimum. 

In Fig. 3(a), we show T '((/>) / "F'(<f> - 0) as a function of <p 
at some T values. There, T' has a minimum at <t> — 0° for 
T > 0.913. That is, the F-phase is stable at high-T and low-// 
near T c , For 0.908 < T/T c < 0.913, the minimum orientation 
<p changes gradually from to 30°, indicating that the L-phase 
is stable. For 0.46 < T/T c < 0.908 at higher fields, the I-phase 
of <p = 30° is stable. The T dependence of the stable orienta- 
tion (f> is presented in Figs. 3(b) and 3(c). At a low temperature 
T/T c < 0.48 at high fields, we see the reentrant transition to 
the L-phase where <p decreases from 30°. These successive 
transitions of the vortex lattice, i.e., L — > I — > L — > F from 
a low-T along H C 2, qualitatively reproduce the phase diagram 
observed in the SANS experiment. 16 Interestingly, we note 
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Fig. 3. (Color online) (a) T'((/>)/T'((f> = 0) as a function of </> at T/T c = 
0.908, 0.909, 0.91 1, 0.912, and 0.913, showing transition I -> L -> F, for (i) 
Vcr,n = 0.15VCT.O-, V niJI = 0.32VD-O-. The intervals of data points are 1° along 
each line. Arrows indicate the minimum of each line, (b) Stable orientation tj> 
of triangular vortex lattice as a function of T/T c along H c 2 line. In addition 
to case (i), we also show lines for (ii) Va- ?7r = Vg-^/3, V KtX = 0, (iii) V a ^ = 
W2, V nJT = 0, and (iv) V a „ = V^/3, V nJT '= V^JL (c) The high-T 
range of the transition I — > L — > F of (i) is focused on. 



that data points at the lowest T (2 K) and high field (1.8 T) 
show (p ~ 23°, 16 which nicely correspond to our result of <p at 
the lowest T in Fig. 3(b). 

To confirm that these qualitative results in Fig. 3(b) do not 
strongly depend on the interaction parameters within appro- 
priate values for MgB2, we show stable orientations also in 
other cases (ii) V^„ = V^/3, V n , n = 0, (iii) V^ = Va-^/2, 
V^ = 0, and (iv) V^ = V^A V„, n = V^/2. The gap ra- 
tios are AJA a =0.24(i), 0.29(h), 0.40(iii), and 0.57(iv). We 
see similar successive transitions there, and find that the range 
of the F-phase (I-phase) becomes wider (narrower) with the 
increase in A^/A,,-. This indicates that the stability of the F- 
phase is due to the contribution of the 7r-band. 

To clarify the contributions of the two bands, which are 
competing with each other, we evaluate the stable orientation 
of each band. As shown in Fig. 4, if we assume that only the 
cr-band is superconducting (Vo- t „ = V n . n - 0), the stable vor- 
tex lattice is the I-phase with <p - 30° at T > 0.45r c , re- 
flecting the Fermi velocity anisotropy in Fig. 2(b). Namely, 
the nearest neighbor of vortices is oriented in the direction 
where D„.Qp v ) is small. At a lower T and a higher H along 
H C 2, (p decreases from 30°, changing to the L-phase. On the 
other hand, if we assume that only the 7r-band is supercon- 
ducting (Vo-.tt = Vtj^a- - 0), we obtain opposite results, re- 
flecting the opposite Fermi velocity anisotropy in Fig. 2(c). 
That is, the stable vortex lattice is the F-phase with </> - 0° 
at T > Q.65T C , which changes to the L-phase with increas- 
ing <p from 0° at a low T. The reason why the vortex lattice 
orientation moves at a low T along H C 2 in both cases is the 
Fermi velocity anisotropy changing its contributions at high 
fields. That is, the anisotropy in the a (7r)-band contributes to 
stabilize the I- (F-) phase at low fields, but makes the low- 
field-phase unstable at high fields. This effect is also seen in 
the case of a square vortex lattice of four-fold symmetric su- 
perconductors. 1 To discuss the reason for these behaviors, we 
consider the yS-model of the cylindrical Fermi surface. : If we 
assume sixfold anisotropy |v,| oc 1 + 0.4cos6^>,, for the Fermi 
velocity in Eq. (2), the stable orientation is at 0°. Instead, if we 
assume anisotropy Di(ep) oc 1/(1 + 0.4cos6<y? v ) for the Fermi 
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Fig. 4. (Color online) Stable orientation <p as a function of T/T c along H^i 
line when only the cr-band is superconducting (Va- M = V„ i7l = 0) and when 
only the 7r-band is superconducting (y„- M = V^a- = 0). 



using the first-principles band calculation to obtain valuable 
information on the multiband and anisotropy of superconduc- 
tivity from the vortex lattice configuration and morphology. 
We hope that this approach will be applied to other supercon- 
ductors in future studies. 

We thank M.R. Eskildsen for stimulating discussions and 
information on their experimental results. We also acknowl- 
edge useful discussions with P. Miranovic, N. Nakai, and 
H.M. Adachi for the formulation and calculation methods. 
This work was supported by KAKENHI No. 21340103. 



surface DOS, the stable orientation is at 30°. Combining both 
anisotropies of v, and Dj(ep), the stable orientation is at 0° at 
low fields, and changes to 30° at high fields. Since there is 
usually the relation D,(£f) cc l/|v,|, the anisotropies of v, and 
Dj(e v ) tend to compete with each other, inducing the change 
of stable orientation between low fields and high fields. 

On the basis of the results in Fig. 4, we discuss the origin 
of the phase diagram of the vortex lattice presented in Fig. 
3(b). The F-phase near T c is given by the contribution of the 
7r-band, since the total DOS N„ and the anisotropy of the n- 
band is large compared with that of the cr-band, as discussed 
in Fig. 2. However, the contribution of the 7r-band decreases 
at high fields, since the superconductivity in the 7r-band be- 
comes normal-state-like there. 7,13 This is because the effec- 
tive upper-critical field of the 7r-band is smaller than that of 
the cr-band, which is roughly proportional to the gap ampli- 
tude Ajt/Aq-. Thus, since the effective upper-critical field de- 
creases for (iv) — > (iii) — > (ii) — > (i), the 7r-band contributions 
becomes weaker in this order in Fig. 3(b). The I-phase ap- 
pears at low-T and high fields, owing to the contributions of 
the cr-band. Therefore, the transition I — > L — > F in the high-T 
region occurs owing to the competition of sixfold anisotropy 
between the two bands. The low-T reentrant behavior from 
the I-phase to the L-phase in Fig. 3(b) is due to the anisotropy 
of the cr-band. The reentrance to the L-phase occurs similarly 
in the vortex lattice behavior of the cr-band in Fig. 4. 

When our results are compared with those in experimental 
observation, 16 it is noted that the phase boundaries in Fig. 1 
are extrapolated near H C 2- The weak anisotropy of the s-wave 
superconducting gap and impurity scattering, which are not 
considered in this study, may be factors of minor quantitative 
modification for our theoretical estimation of a stable vortex 
lattice configuration. 

A recent SANS experiment 29 on CaAISi, which has the 
same Allotype crystal structure and a similar two-band fea- 
ture, 30 has shown a first-order orientation transition of the vor- 
tex lattice. The present methodology is readily applicable to 
this and other omnipresent multiband superconductors, such 
as Fe pnictides. 

In summary, a method for the microscopic Eilenberger 
analysis of a stable vortex lattice configuration along H& 
combined with the first-principles band calculation of DFT is 
developed for a two-band superconductor. We apply this new 
formulation to MgB2, and explain the successive transition 16 
L— > I — > L — > F in the orientation of a triangular vortex lat- 
tice, by properly considering the Fermi surface anisotropy of 
the two bands. This successful endeavor establishes a method 
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